Study on two coupled relativistic superfluids with spontaneous 

symmetry breaking 



Sun Zhang 

Department of Physics, Nanjing University, Nanjing 210093, China 

Abstract 

In this Letter, we have studied two coupled relativistic superfluids with spon- 
taneous U(l) symmetry breaking, using Poisson bracket technique. After con- 
structing the commutators between thermo-quantities and field quantities, 
the equations of motion are obtained. These equations describe the system 
in the frame of two-constituent superfluid theory and provide a clear picture 
relating the symmetry breaking and the physical quantities explicitly. The 
interference effect is discussed and the Josephson-type equation is also given. 
Furthermore, the dissipation effect and entropy production mechanism are 
presented and the dissipative coefficients are given explicitly in the first-order 

theory framework. 

PACS: 47.75.+f; 47.37.+q; 67.40.Bz 
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In 1941, Landau proposed his famous two-constituent theory on superfluid 4 He |],@. The 
central idea of the theory is that at temperature below the critical point, helium behaves as 
if it were a mixture of two different liquids. One of these is a superfluid, and moves with zero 
viscosity along a solid surface, while the other is a normal viscous liquid. No momentum is 
transferred from one to the other. The two-constituent superfluid theory was studied and 
generalized by Khalatnikov || brilliantly in the non-relativistic frame. 

The first approach to the relativistic superfluid theory was provided by Israel |§] and 
Dixon m concerning the perfect fluids based on the idea of two-constituents originated 
by Landau. The medium is usually called relativistic in two senses: if it has a relativistic 
equation of state or when it flows at a relativistic velocity. While these conditions can be 
both satisfied in the neutron stars whose interior construction is composed of a superfluid 
nuclear matter ]7| . The method proposed by Israel and Dixon is useful for some calculations, 
but arising a problem. 

Generally speaking, the coupled two-constituents are not perfect fluids strictly and the 
superfluid component cannot be regarded as a completely independent fluid H . That is, the 
coupling effects lead the deviation from ideality ||. It is meaningful to take into account the 
deviation from the perfect fluids. Khalatnikov and Lebedev || and Carter and Khalatnikov 
Tcfl suggested two equivalent approaches to include the interaction between the superfluid 



and normal fluid. 

Furthermore, Son [[U]] presented another approach to the relativistic superfluid with 



U(l) symmetry breaking based on the Poisson bracket technique ||12 |. Such a description 



is also arisen from the two-constituent theory, but it is more clearer to show the relation 
between superfluid and the symmetry breaking. After constructing some basic relations of 
the fundamental quantities by Poisson brackets, the dynamic equations are obtained directly, 
which provide the explicit meaning about the symmetry breaking. 



It is possible to generalize the discussion presented in Ref. |11|] in several directions. The 
most important one is that we should construct a clear frame to discuss the effects between 
two coupled relativistic superfluids including the coherence effect and the dissipation effect. 



For the symmetry spontaneous breaking and the phase coherence are assumed from the very 
beginning, such an approach could be looked as a natural framework for investigating the 
inter-relationship of two coupled relativistic superfluids in the hydrodynamic limit. 

To study the two coupled relativistic superfluids 1 and 2 with U(l) symmetry breaking, 
the order parameter is of the form (ipi) = \(tpi)\ e"^ (i = 1, 2), the density of the U(l) charge 
Pi = —i(i[)*doi/ji — doi/j*i[>i) and the phase 0« are introduced IflTf] . 

In order to establish the dynamic equations, the Poisson brackets should be written 
down to give the commutative relations between dynamic variables. There are two kinds of 
Poisson brackets: the fluid (thermo-) ones and the field ones. 

The fluid ones are Oil 



[T 0l (x),A(y)] = A{x)di6(x - y), A = p u p 2 ,s (1) 
[T 0t (x), B(y)} = -d, t B(x)5(x - y), B = fa, 2 (2) 
[T°*(x),T°%)] = ^(x)A_ T oi (y) ^_^ (x _ y); (3) 

where T° l is the momentum density in the 4-dimensional energy-momentum tensor, s the 
entropy density. pi ;2 and ipi^ are just the density of charge and the phase of the order 
parameter respectively. 

The Poisson bracket for the field part includes the commutator of pi with 0j (up to an 
adjustable factor) 

[Pi(x),<A?(y)] = -6{x-y)6ij. i,j = 1,2 (4) 

To get the equation of motion, we should construct the most general Hamiltonian de- 
pendent on the fluid quantities and the field quantities 

H = J drT 00 (s,p 1 ,p 2 ,T 0i ,d i( j )l ,d i( j )2 ), (5) 



where T 00 is a functional of dicfii^ but not <pi^ itself due to the U(l) invariance |p~3| , pT4 
The form of the T 00 cannot be identified at this moment, it could be calculated from the 
microscopic theory. 



Then the equation of state, describing the relation of the pressure p, the temperature T, 
the chemical potentials /ii j2 , and the phases <p 1)2 , could be written in the differential form as 

dp = sdT + p 1 dfx 1 + p 2 dfi 2 + V&d Q(^0i) 2 J + V 2 2 2 d Q(a M 2 ) 2 J - (V 2 + V 2 \)d Q^i^&J 

(6) 

where, and V^ 2 2 are the factors characterizing the strength of the internal energy for 
superffuids 1 and 2 respectively, while V± 2 and V 2 \ are the coupling constants between two 
superfluids due to interference and V± 2 = V 21 . In Eq. (|j), the first three terms are fluid 
terms, they describe the normal fluids part, while the gradients of the phases in the left 
terms are specific to the superfluids and have the significance of the superfluid velocities. 

Based on the Poisson brackets (^) — (Q) and the thermodynamic relations, we obtain the 
dynamic equations 

d, (u» Pl - vAdy, + l -{v? 2 + v 2 \)d^ = o, (7) 
d, Lrpi - v 2 \d^ 2 + ~(v* + v 2 \)d^ = o, (8) 

vPdrfa + in = 0, (9) 



« M ^ 2 + /i 2 = 0. (10) 

It is clear that Eqs. (|7]) and (|8|) are the conservation laws of the U(l) charge in the Lorentz 
covariance form. Within the bracket of Eq. (0) [or Eq. (|8|)], the first term presents the 
normal current and the 4-dimensional velocity is the velocity of the normal component; 
while the second term gives the superfluid current and the gradient of the phase is clearly 
the velocity of the superfluid component, and the two-constituent description is constructed 
explicitly. The third term is the interference current, which illuminates the phase coherent 
influence of the superfluids from each other. Then the U(l) charge is a sum of the normal 
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current, the superfTuid current, and the interference current due to the coupling to the other 
superfluid, and Eqs. (|7|) and give the conservation relations of the U(l) charge. 

On the other hand, Eq. (g) [or Eq. (|i~0|)1 gives the time dependence of the condensate 
phase on the chemical potential. The chemical potential in these equations is just the 
change in energy density per change in the charge density, and also the minimum energy 
required to add a paticle to the system. The chemical potential also has the meaning of the 
background field due to the U(l) gauge transformation, and coincides the scalar potential 
V in the superconductivity. If we suppose that the two superfluids are kept at uniform 
chemical potentials, it is direct to obtain the time dependent relation of the difference of the 
condensate phases according to Eqs. @ and fllCD 

u"0 M (0i - fa) + (/ii - /x 2 ) = 0. (11) 



It is just the relativistic version of the Josephson equation, in BEC coherence effect [|15| and 
superconductivity |16[ . 



We can also construct the energy-momentum tensor as a sum of the thermo (fluid) 
part and field (superfluid) part as follows 

= (e + p)v?u u - pg^ + VftPfadTfa + Vi 2 dy 2 d u cj) 2 - V&drfaVfa - V 2 \d^fad v fa, 

(12) 

which satisfies the conservation law 

d,T^ = 0, (13) 



where the fluid part 



TZ d = (c + P)u»u» - pg^ , (14) 



describes the normal relativistic fluid outside the superfluid condensation phase and the 
metric tensor is 

<T = diag(l, -1,-1,-1), (15) 
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in accordance with Landau and Lifshitz, and u^u^ = 1. While the field part describes the 
coherent motion and the interference effect of two superfluids. 

Next, we would like to explore the dissipation effect and the entropy production mech- 
anism of two coupled relativistic superfluids system to provide a full understanding of the 
behaviors described here, following Landau and Lifshitz [0 and Weinberg JT7 . 



Generally speaking, the dissipative processes are important, when space-time gradients 
of hydrodynamic quantities in the system become large relative to its relaxation scales. 
We would analyze the dynamics of small departures of the system from equilibrium states 
and manage to give the dissipative equations as linear responses to the deviations from 



equilibrium states and two famous dissipative theories by Eckart [18| and by Landau and 
Lifshitz would be the special cases in the theory presented here. 

When we study the imperfect system of two coupled superfluids, the dissipative effects 
should be taken into account by modifying the energy-momentum tensor T^ u by a term 
AT^ V , and could be written in matrix form as 



p)u^u u — pg 



hi/ 



-V 2 



21 



V 2 



A ( 

12 
2 

22 



d 1 - 
d 1 ' 



+ AT"", 



(16) 



satisfying the conservation equation as 



dJT^ + AT"") = 0. 



(17) 



The other dissipative effects for the charge and the phase could be considered as 



and 



0„ 



u fl p 1 + Apt 



\ 



u^p 2 + Ap£ 



d„ 



( A \ 

<Pi 



( 

V 



v 2 



-V 2 

V 12 



-V 2 V 2 

v 21 v 22 



(18) 



Pi + Acpi 
,/i 2 + A(f) 2 j 



(19) 



The dissipative terms should satisfy the following constraints 

u.u v AT^ = 0, 



(20) 



(21) 



For comparison, the constraint used by Landau and Lifshitz is u u AT^ u = and that of 
Eckart is Aff = 0. 

Then using the equations of motion and the thermodynamic relations, we can get the 
entropy production equation arisen by the dissipative terms, 



/ A ,\ 



y 



y2 _y2 \ I o„ , \ 



-T/ 2 T/2 
"21 v 22 



d^(j) 



\ d^2 y j 



+-AT^u v - —^AT^dvT. 



Where the entropy density is defined as 



Apr 

^ y 



By means of the projection operator (g^ — u^u"), AT' 11 ' can be expressed as 

AT» V = {{g^ a - u»u a )u v + (g va - u v u a )u») q a 

+ {g» a - u»u a ){g u P - u v u^)T aP + (g^ - u»u v )t. 

Defining the quantities as 



/ AA, \ /_ _ \ ( „\ 

T 



AN i 
AN 2 



\ 



J 



V T J 



(22) 



(23) 



(24) 



(25) 



0\\ G\2 
021 &22 

where o\\ and 022 are the charge diffusion coefficients for superfluids 1 and 2 respectively, 
while a 12 = &21 are the mutual diffusion coefficients between two superfluids and are sup- 
posed to have only nonnegative eigenvalues. Then the charge dissipative equations are 



Aft 



{<f - u»u u )d u 



<?\\ &12 



C21 0"22 



LV 

\ AN 2 J 



T 



(26) 



Then the first term in Eq. could be written as 



U t* [ rp 5 y 



Apt 



■or > f)^ 



yAN 2 j 
( 



-W-u^d^AN! , AN 2 ) 



\ 



\ 



1 < m \ 



J 



\AN 2 j 



^11 012 
0"21 °22 

where the positivity of the inverse matrix of a is assumed. Defining the quantities A± 2 as 

\ 



(27) 



V 2 -V 2 

v ll v 12 



-V 2 V 2 

v 21 v 22 



and AMf 2 as 



^ AMf ^ 



y AM 2 j 



( 



J 



\ 



(28) 



\ A 'l 



(29) 



C11 C12 

C21 C22 

where C11 an d C22 are the bulk viscosity coefficients for superfluids 1 and 2 respectively, 
while C12 = C21 are the mutual viscosity coefficients between two superfluids and have only 
nonnegative eigenvalues supposedly 

Then the second term in Eq. ( P2"D could be written as 

\ 



A0i A02 

T ' T 

A(j) X A(f) 2 

T ' T 



-V 2 V 2 

v 21 v 22 



V 



V 2 -V 2 

v ll v 12 



-d, (AM 1 , AM 2 ) 



Cn C12 

C21 C22 



\ 



1 / \ 

' AM" x 



0„ 



y AM 2 j 



(30) 



if we present the dissipative equation A0 ii2 as 



A0i 


( c 


c ^ 
















v A< ^ 2 J 


^ C21 


C22 J 









(31) 



where the viscosity coefficients Ci an d C2 axe those related to the normal velocity of the 
systems respectively. The dissipative components in AT^ of Eq. (|24| ) could be shown as 

1 



kT 



T 



(32) 



Ta/3 = V(d a Uf3 + 3pU a - -gapd^), 



T 



(33) 



(34) 



where a is the charge- diffusion coefficient, £ is the total bulk viscosity coefficient, while k 
and 77 are the thermal conductivity and the shear viscosity respectively. According to the 
Onsager's principle, 



& = & i = l,2 

and the positivity should be ensured by the following inequality 

C(CllC22 — C12C21) ^ C2C1I — 2ClC2Cl2 + C1C22, 



(35) 



(36) 



such an inequality is a generalization of that by Landau and Lifshitz || . 

Then the full set of dissipative equations are established here, and all transport coeffi- 
cients are defined. These coefficients could be calculated from the microscopic theory and 
Kubo's linear response theory could be used here to build up the connections between the 
microscopic fluctuation and the macroscopic dissipation. 

In summary, the fundamental equations of motion have been constructed which could 
be used to study the effects between two coupled relativistic superffuids. The theoretical 
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frame presented here are essentially based on the idea of two-constituent model proposed by 
Landau. Furthermore, in such a theoretical description, there arises more clearer physical 
significance for the equations of motion, which relate the symmetry breaking and the physical 
quantities naturally Based on these, the interference effect and the Josephson-type equation 
are discussed; the dissipation effect and entropy production mechanism are presented, and 
the dissipative coefficients are given explicitly in the first-order theory framework. 

Besides the interference effect between two coupled relativistic superfluids with sponta- 
neous U(l) symmetry breaking, another important effect is the discontinuity, including the 
shock wave propagation in the two-constituent relativistic superfluids. Although the shock 
wave was studied by many authors in the relativistic systems [pl|, a new theoretical frame 
is still necessary to provide some new features qualitatively. Instead of working in the frame 
of perfect fluid, we would rather select the theoretical picture presented here to give a new 
description on the shock wave in relativistic superfluids [EO]. 
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